A comprehensive study of electric, thermoelectric and thermal conductivities of 
Graphene with short range unitary and charged impurities 
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Motivated by tlie experimental measurement of electrical and hall conductivity, thermopower and 
Nernst effect, we calculate the longitudinal and transverse electrical and heat transport in graphene 
in the presence of unitary scatterers as well as charged impurities. The temperature and carrier 
density dependence in this system displays a number of anomalous features that arise due to the 
relativistic nature of the low energy fermionic degrees of freedom. We derive the properties in detail 
including the effect of unitary and charged impurities self-consistently, and present tables giving 
the analytic expressions for all the transport properties in the limit of small and large temperature 
compared to the chemical potential and the scattering rates. We compare our results with the 
available experimental data. While the qualitative variations with temperature and density of 
carriers or chemical potential of all transport properties can be reproduced, we find that a given set 
of parameters of the impurities cannot fit all the observed data. 



I. INTRODUCTION 

An unusual new electronic structure and the possibil- 
ity of graphene as the basis for technologies of the future 
has sparked an intense effort in its fabrication and char- 
acterization. A number of spectacular properties such as 
conductivity in the limit of zero carrier dcnsitji"— , per- 
fect tunneling through potential barriers^ri and quantum 
hall effect at room temperatures^!^ have already been ob- 
served. While the anomalous properties of electrical con- 
ductivity has received a lot of attention, data on ther- 
mopower, Hall conductivity and Nernsli^"— have now 
revealed temperature and gate voltage dependence which 
need to be understood consistently within a single trans- 
port theory. Excellent summary of the previous theoreti- 
cal and experimental work has recently become available 
while this work was in progres s ^'^i^^ . 

Graphene is a two dimensional allotrope of carbon with 
a hexagonal crystal structure. Since it is made up of two 
interpenetrating triangular sublattices, the unit cell has 
two atomsi^. As long as the sublattice symmetry is pre- 
served the two bands touch at two points in the Brillouin 
zone. In the vicinity of these points the hamiltonian is 
linear in momentum and has the structure of k • ct where 
k is the momentum and a = {a^, cr^, a'^]^. The Pauli 
matrices represent pseudospin with the two components 
referring to the two sublattices. The linear dispersion 
means that the electrons near these points in the Bril- 
louin zone behave like relativistic masslcss particles in 
the absence of impurities and arc called Dirac points^i^. 
For pure graphene the fcrmi surface is at the Dirac point 
and the density of states depends linearly on energy near 
the chemical potential. The vanishing density of states, 
the conservation of the operator k • a and the existence 
of two zeroes (or equivalently) valleys in the band struc- 
ture are responsible for a number of novel phenomena in 
graphene. 

The most striking observation is that the electrical con- 
ductivity varies linearly with carrier density when the 
carrier density is not too small and that it is nonzero 



even when the carrier density goes to zeroi"— . Numerous 
attempts at explaining the latter has led to a number of 
different values for the minimum conductivit y ^^' ^^i ^^~ '—. 
The reason for the different predictions can be traced to 
sensitivity of the results to the different approximation 
schemes and order of limits employed in the calculations. 
For example taking the zero frequency limit before the 
zero temperature limit docs not commute with the limits 
taken in the opposite order. Introducing an additional 
scale, such as the scattering rate further complicates the 
order of limits providing a wide spread of possible values. 
Experimentally it is clear that the observed minima is 
not universal. Within a Boltzmann transport formalism 
a scattering rate inversely proportional to the energy can 
account for the observed linear dependence with respect 
to carrier density. One possible source of such scattering 
is long range Coulomb scatterersi^ii^i^ii^. While this 
theory works well for finite densities, the finite minimum 
conductivity requires new physics near the Dirac point. 
Based on the observation of charge inhomogeneity in this 
limit^"— , a possible resolution is that the Coulomb scat- 
terers promote the formation of charge puddles. These 
puddles mask the approach to zero carrier density and 
provide an effective mechanism for minimum conductiv- 
ity. Another possibility which we work out is due to the 
fact that for small charge densities, the effect of Coulomb 
scattering due to point charged defects also needs to be 
calculated self-consistently. 

Alternatively, a mechanism that can provide a sim- 
ilar dependence of scattering rate on energy is strong 
scatterers in the unitarity limit^i^. Within this ap- 
proach, the scatterers introduce resonances and, in the 
independent scattering approximation, an effective impu- 
rity band forms which provides a finite density of states 
in the vicinity of the node. The width of this impurity 
band is set by the density of scatterers and for energies 
larger that the impurity bandwidth the linear density of 
states is recovered. Crucially the same parameter, i.e. 
density of impurities sets both the band width and the 
scattering rate. Qualitatively a constant conductivity at 



low densities crossing over to a linear in carrier density 
behavior is expected. A similar result occurs for Coulomb 
scatters as well, since for low carrier densities, even weak 
impurity potentials can induce resonances at low energies 
as the density of states goes to zero. 

In this paper, we provide the dependence of longitu- 
dinal and hall conductivity, thermopower, thermal con- 
ductivity and Nernst coefficient for various regimes in 
temperature, scattering rate and chemical potential for 
unitary scatterers as well as calculate the self-consistent 
scattering rate for Coulomb scatterers so as to compare 
the two cases. For a finite impurity bandwidth, the lead- 
ing contribution to conductivity at the node is indeed 
universal in the limit of zero temperature. For temper- 
atures smaller than the impurity bandwidth, the correc- 
tion to the universal value is of order (Tr/Zi)^, where r is 
the mean free path. We compare our results to the only 
data available where conductivity, Hall resistance, ther- 
mopower and Nernst were measured on the same sample. 
Our chief conclusion is that no single choice of impurity 
concentration can account for all the observed data. For 
unitary impurities, none of the data can be fit in the en- 
tire range of gate voltages measured. For charged scat- 
terers, the longitudinal conductivity can be fit but not 
the thermopower or Nernst. Since the experimental data 
is obtained from a two probe measurement, the fact that 
the Hall resistance cannot be reproduces is not surpris- 
ing. 

In section II, we derive the form of the impurity av- 
erage self energy and discuss the nature of the effective 
dispersion and mean free path. A similar analysis for long 
range Coulomb impurities is presented in section HI. The 
transport formalism used to derive the conductivities is 
discussed in section IV. The form of electrical conductiv- 
ity and hall resistance is analyzed in section V. Section 
VI discusses the nature of thermoelectric properties. In 
section VII, we present the results for thermal transport. 
We compare our results for conductivity, hall coefficient, 
thermopower and Nernst with experimental data in sec- 
tion VIII. The details of all calculations are available in 
the appendix. Many parts of our work have already ap- 
peared in separate works of many others as we note in 
the References; our contribution is primarily the compre- 
hensive calculations and comparison with experiments of 
diverse transport properties. We present asymptotic an- 
alytic expressions for the various transport quantities in 
a set of tables; these may be especially useful since they 
readily provide a physical basis for the experimental re- 
sults. 



II. IMPURITIES SELF AVERAGING 
FORMALISM AND CARRIER DENSITIES 



The Hamiltonian for graphene in tight binding formal- 
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Where t is the nearest neighbor hopping amplitude and is 
related to the fermi velocity by Wi? = | ^ . The operators 

in the Hamiltonian {cLl,b'-} represent electron creation 
operators on sites i and j in the graphenc's honeycomb 
lattice which belong to the A and B sublatticc respec- 
tively. The two atoms per unit cell leads to a 2 x 2 
matrix for the Green's function of graphene. The two 
bands touch at two points in the Brillouin zone labelled 
by K and K . In the vicinity of these points the Greens 
function is^ 
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Where the function Q{ k ) is equal to —^ + arg{kx -t- 
iky) and the dispersion relation at the node is given by 

(j){k) = ±hvF\ A: |. In the presence of impurities we have 
an additional term in the Hamiltonian given by 
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In this section, the model chosen is s-wave scatter po- 
tentials in the unitary limit. If the assumption is made 
that we have identical impurities randomly distributed 
throughout the system, an impurity self-average Green's 
function can be used as an approximate solution to the 
full Green's function of graphene. The impurity self- 
averaging is valid if the sample size is much larger than 
the coherence length of electrons which is the case for 
most experiment. In the dilute impurity limit, the scat- 
tering of single impurities dominates transport. The self 
energy is calculated in the full self-consistent Born ap- 
proximation (FSB A). It is important to emphasize the 
regime of validity of thus approximation. The FSBA 
is known to fail due to interference effects that involve 
scattering of multiple impurities^i^ when the chemical 
potential is at the Dirac point. However, at finite chem- 
ical potentials the approximation is justified as long as 
weak localization effects are not important^. Further 
support for the FSBA approach is seen in the agreement 
of the density of states determined within this approx- 
imation and exact numerical methods [see |3J and ref- 
erences therein]. We use the FSBA to study transport, 
our results are valid for finite chemical potentials but will 
ultimately breakdown when the chemical potential is at 
or near the Dirac point. Analysis of where the crossover 
occurs is beyond the scope of this study. 
The effect of the impurity states on transport quantities 
is to produce a finite value for conductivity at the node as 



the carrier density is no longer zero. In the dilute limit of 
impurities, the self energy in the full self-consistent born 
approximation is given by22 
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For s-wave short range scatterers, the self energy is mo- 
mentum independent. For Coulomb scatterers, we will 
assume momentum independence of the self energy. This 
assumption breaks down near the Dirac pointSl. The 
electronic Carrier density and the change in carrier den- 
sity near the node due to a change in chemical potential 
/x are given by 
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where {Ny,Ns} represent graphene's valley and spin de- 
generacies and G^ are the retarded Greens function. The 
integral over momentum can be performed leading to 
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where D is the electronic band width and the functions 
{A{e),B{£)} are given by {e - i?eE(£), -ImT.{e)). 



A. Self Energy 

In fig[l] we plot |/toI]| as a function of energy (e) for 
several impurity concentrations. The impurity concen- 
tration is defined as the number of impurities per Car- 
bon atom per unit area. As the impurity concentra- 
tion is increased the scattering rate (which is propor- 
tional to |J?7iE|) increases. It is weakly dependent on 
energy near the node crossing over to a 1/e dependence 
for large energies. The crossover scale is determined by 



|ImZ| 
0.08 p 



006-1 



' / 0.04 






/ 



,,''' / ''/ °^ 



"Imp 




6.7 X 


10-' 





1.2x 


10-* 





1.7x 


10-* 





4.1 X 


10-* 




6.5 X 


10-* 









FIG. 1: The absolute value of the imaginary part of the self 
energy is plotted as a function of the energy. The impurity 
concentration of unitary scatterers used to find the self energy 
curves is given in the legend. 



the bandwidth (D) times the square root of the impu- 
rity concentration^. This crossover scale is the impurity 
band width. One can regard the effect of unitary scat- 
ters as producing resonances that fill in the density of 
states in the vicinity of the node. The linear density of 
states is recovered beyond the impurity band width. Un- 
like Born-scattering, beyond the impurity band width the 
scattering rate remains inversely proportional to the en- 
ergy. Qualitatively, in this regime, the physics is similar 
to having weak Coulomb scatterers. Physical insight is 
gained by taking various limit of the change in electronic 
carrier density with respect to temperature, chemical po- 
tential and impurity band width. 
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In the limits above, the electronic band size has been 
taken to be the largest energy scale and the real part of 
the self energy absorbed in a suitable redefinition of the 
chemical potential. In obtaining eqnllH we have assumed 
that the imaginary part of the self energy is roughly con- 
stant up to the impurity band width with its value deter- 
mined at zero energy. The integrals are approximated as 

/&^(...)^r;:,#(---)and/de^(...)« 
/ de6(p — e){- ■ • ) = (• ■ • )|£=^ at high and low tempera- 
tures respectively. Since 6n = J d^li. J^ dedkS {e — e(k)) 
we can extract an effective dispersion relation in the dif- 



ferent regimes 
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In the limit where \fj,\ >> \Im'S\ >> T, the Dirac disper- 
sion relation is preserved. In the other limits, the disper- 
sion is effectively that of a free electron with a mass de- 
termined by the dominant energy scale. Given this form 
we can use the effective dispersion in eqnlTHto calculate 
the longitudinal conductivities in various regimes^^. 
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For small chemical potential, we find that the conductiv- 
ity is constant and depends logarithmically on the imagi- 
nary part of the self energy which in turn is proportional 
to the square root of the impurity concentration^^. Thus, 
within this picture, we obtain a finite minimum conduc- 
tivity which is not universal. As we increase the chemical 
potential, at low temperatures, we crossover to a regime 
where the conductivity scales as /x^ oc 6n (note in this 
regime /toE cx l//i). Both these features (non- universal 
minimum conductivity and linear dependence on carrier 
density) are in qualitative agreement with observed data 
in graphene. There is also a constraint implicit here that 
the crossover scale is controlled by the same parameter 
that determines the value of the conductivity minimum. 



B. Mean Free Path 

We study the transport properties in weak magnetic 
fields where we are in the hydrodynamic limit (uJcT ^1). 
In figl2]we plot the ratio of the mean free path to the cy- 
clotron radius as a function of the chemical potential for a 
fixed magnetic field. The mean free path is roughly con- 
stant up to the impurity bandwidth which in this case 
is 0.06 eV. Beyond this scale, the mean free path grows 



suggesting an energy dependent scattering rate that be- 
comes smaller at higher energies. This behavior is con- 
sistent with the behavior of the imaginary part of the self 
energy. For the majority of the range shown in figlU the 
mean free path is much smaller than the magnetic length 
for a field of 8 Tesla, for this range we are within the 
hydrodynamic regime. 
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FIG. 2: The mean free scattering length (Limp) in units of 
cyclotron radius (Lb) is plotted as a function of the chemical 
potential (/i). The impurity concentration of unitary scat- 
terers is rii = 2.2 x 10"* and the cyclotron radius is calcu- 
lated in an 8 Tesla magnetic field. The Fermi velocity used is 
vf — 1.0 * lO^m/sec 



C. Specific heat 

Thermoelectric transport coefficients, such as thermo- 
electric power, Nernst and thermal conductivity depend 
on the specific heat at constant volume. The energy den- 
sity and specific heat in graphene in the presence of uni- 
tary scatterers are given by^ 
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With the same approximations used to derive the con- 
ductivity above, the specific heat and electronic carrier 
density, to leading order in the appropriate small param- 



eter, in different regimes are given in tablelH The specific 
heat divided by temperature is proportional to dSn/dfj,. 
This is natural as we expect the two quantities to be 



TABLE I: Carrier Density and Specific Heat 



Quantity 
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proportional to tiie density of states which in turn is the 
imaginary part of the self energy. 



III. COULOMB SCATTERERS 

In this section, we consider the nature of scattering 
in the presence of charge impurities. We follow the same 
approach as in the case of unitary scatterers and compute 
the self energy in the self consistent Born approximation. 
To first order, within born approximation the self-energy 
due to screened Coulomb scatters has the form: 
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where rii is the concentration of charge impurities, e is 
the charge of an electron, n is the area, eg is the vacuum 
permittivity, G^{k' ,iuin) is the green's function in the 
absence of impurities and k is the permittivity of the sub- 
strate. The function qTpi k , iw„) is the inverse Thomas 
Fermi screening length. In general, the self energy is 
a 2x2 matrix and must be handled in a self-consistent 
manner. In order to correctly account for changes in the 
ground state energy the Green function in ea. (|18p must 
be replaced by the full Green's function. 

For k = the self energy is diagonal and has the form: 
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The density of states and the inverse Thomas Fermi 
screening length are given by 



To obtain the transport coefficients in the presence of 
Coulomb scatterers we substitute the scattering rate ob- 
tained here in expressions derived for unitary scatter- 
ers. An important point to emphasize is that the self 
consistency yields a finite screening length at the node 
while crossing over to the inverse fermi wave vector at 
large carrier densities. Thus, the Coulomb potential 
is screened providing a mechanism for finite density of 
states and conductivity at zero bias. An alternative ap- 
proach to obtain effective screening is to build on the 
observation of inhomogcneities near the Dirac point^^. 
The real space density fiuctuations lead to screening 
that can be determined by solving for the Thomas-Fermi 
screening length self consistently given a random poten- 
tial distribution22i2S. Determining whether the impurity 
states obtained via our self consistent approach is equiv- 
alent to the effective inhomogencous electron gas in real 
space requires is beyond the scope of this article. 

The self energy is momentum independent at the node. 
We assume that vertex corrections are small and compute 
the self energy considering the mass-shell only. While 
this is exact for unitary scatterers and for large chemical 
potentials for Coulomb scatterers, the nature of momen- 
tum dependence of the vertex correction near the node 
for the latter case is unclear. While a sizable momentum 
dependence has been calculated^, the study ignores the 
effect of screening due to the formation of mid gap states. 
Our approach is to solve for the effective screening deter- 
mined by the self energy which renders the approxima- 
tion of momentum independence applicable. For a finite 
impurity concentration ca. (|20[) is simplified to the form: 
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FIG. 3: The absolute vale of the imaginary part of the self 
energy is plotted as a function of the energy. The impurity 
concentration of Coulomb scatterers used to find the self en- 
ergy curves is given in the legend. 



Qualitatively the imaginary part of the self energy is 
similar to that of unitary scatterers. At large carrier den- 
sities, the scattering rate falls off as Xj [i. This is clear in 
figl3]for smaller impurity concentrations. The divergence 
at zero chemical potential is cutoff by the emergence of 
a finite screening length. The finite scattering rate is re- 
sponsible for the observed minima in conductivity within 
this scenario. The existence of charge inhomogeneities at 
low bias has been experimentally observed and does lead 
finite conductivity at the Dirac point. Whether the real 
space realization of this phenomena yields charge pud- 
dles is an open question that is beyond the scope of this 
work. 

IV. TRANSPORT FORMALISM 



In response to applied electromagnetic fields and ther- 
mal gradients, the electrical and heat current are induced. 
Within linear response formalism, these quantities are re- 
lated as^^i^ 
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The summation over momenta can be performed in 
closed form, but the solution can only be obtained nu- 
merically. The screening obtained at zero energy allows 
us to approximate the self energy to be independent of 
momentum and compute its dependence on the chemical 
potential. The results are shown in figEl We have cho- 
sen impurity concentrations that yield the best fit to the 
data analyzed. The only free parameter is the distance of 
the impurities from the graphene sheet. Unlike unitary 
scatterers, which are part of the graphene layer itself, the 
Coulomb scatterers are in the substrate. In these calcu- 
lations, we have assumed that the charged impurities are 
on the graphene sheet, thus requiring a very small con- 
centration to fit the data. 
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where J is the charge current density, Jq is the heat 
current density, a'^^^ is the electrical conductivity, /S'^' ~ 
/3(2) is the thermoelectric conductivity and k!^'^^ is the 
heat conductivity. The heat current density is related to 

the charge current density : Jq = Je — fJ-J , where Je 
is the energy current density. Each of these conductivity 
tensors is computed using retarded current-current corre- 
lation function within the standard Kubo formalism. The 
current-current correlation function and current densities 
are given by^"— : 
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where ja,i(T,q) is the current density operator, a refers je^aifijt) and JE.aifiTt) are the electrical and energy cur- 
to the type of current density, i refers to the components, rent densities^, {■0, V'^} are the fermion annihilations 



and creation operators and Ar is the vector potential. 
The conductivities can be related to the appropriate cur- 
rent current correlation function. For closed boundary 
conditions, the thermoelectric transport coefficients 5, 
are related to the conductivities /^^^^ and ct^°): 
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where S'a;^: is the thermopower and Cy is the Nerns t^^i^^ . 
The Ncrnst coefficient v is defined similarly to the hall 
coefficient, but here the important quantity is the ratio of 
transverse electric field to longitudinal temperature gra- 
dient. Similarly, the thermal transport coefficients can 
be obtained fro m^^i^^ 



K 



^(3) ^(2)(^(0))-l^(l) 

T T 



(27) 



The components can be related to Nernst and ther- 
mopower: 
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In metals the thermoelectric transport coefficients such 
as Ncrnst and thermopower are not large and generally 
do not contribute to thermal conductivity. The situa- 
tion is more interesting in materials where the density 
of state vanishes, such as graphene and high tempera- 
ture superconductors. The focus of this paper is on the 
anomalous dependence on gate voltage and temperature 
of thermopower and Nernst, the latter being orders of 
magnitude larger than typical metals. 



ELECTRICAL CONDUCTIVITY 



In this section, we discuss the formalism for calculating 
the electrical conductivity tensor and a comparison of 
our results with other theoretical studies is given. The 
electrical conductivity tensor is calculated using Kubo 
formula^S: 
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(29) 



where Ilf 'f is the current-current correlation function and 

. '-^ 
the indices e referring to the vertex corresponding to 

charge current. Using the definition of the charge current 

density operator in eqni25l the resulting current-current 

correlation function is: 
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The calculations arc presented in the appendix. The con- 
ductivity in the presence of unitary scatterers 1 023,35,42 
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Electrical transport 


coefficients 








Quantity 


iI)T, M < 


h 

T 




(//)T, 


h 

T 


< ^x 


(///) 


ft 

r 


< /x <T 


dn 


'^l^-^|i-(|i^sp 


) 

-H 


sgn(/x)|/x|^ 1 


Tn 1 


n^h^vl 


27rft2^^p ' 




2+|(^f + |(^f + --- 


TTif +... 


2 ft ' 4 'J'-^ "•" 


(^xy /CTO 


3 h "T 


sgn(n)-K(wcT)^ , 


2 T "•" 


tan Oh 


16(k;ct)'^ /j,r , 
3 h ~^ 


2 ^ttr^^ 


Krf(^)f + ••• 


Rh 


8{WcT)^ /XT 1 


-sgnifi) {wct)^ . 


'M2 + ... 

^ /J,T / ' 


-1 2(lUeT)^ 


\TtJ t ^ 


3|B|cro h 1 




\B\cro 2tt ' 


|B|cro TT 



V 












^ 


Temperature 


50 




50K 

255 K 


40 




400 a: 



,f 



/ 






The kernels of the electrical conductivity tensor have 
been defined in terms of the functions {A, B} ~ {e — 
Re'E,{e), — /mS(e)} (for details sec appendix A). Defining 

> , wc analyze the dependence 



O-Q 



^ and w;2 



of conductivity in the four regimes described earlier. 



Longitudinal and transverse electrical 
conductivities 



FIG. 4: The diagonal electrical conductivity is plotted in units 

2 

of ^ as a function of chemical potential (/i). The impurity 
concentration of unitary scatterers is rii — 1.7 x 10^'*, while 
the temperature of each curve is given in the legend. 



The dependence of the longitudinal conductivity and 
hall resistance are shown in figH] and figEl Longitudi- 
nal conductivity has a minima at the node crossing over 
to a ^^ (linear in charge density) dependence for large 
carrier densities. The crossover occurs at the impurity 
bandwidth. The hall coefficient is linear at the node and 
falls of as l/fJ,^ for large chemical potential. 



Analytic expressions obtained in various asymptotic lim- 
its are shown in table HIl In the regime where {T << 
1^1 << ^}, the longitudinal electrical conductivity is a 
linear function of fi/\Im'E\. The slope in this regime is 
equal to half the conductivity quanta and agrees with the 
other theoretical results^^'^'^. For unitary scatterers in 
this regime t(^) ^ ^ which implies that the conductivity 
is proportional to carrier density. The Hall coefficient is 
inversely proportionally to (///|/77iS]|)^. The coefficient of 
proportionality depends only on physical constants and 
the scattering rate. Futhermore, the Hall coefficient in 
terms of the carrier density is Rh = — 1/ecn^ agrees with 
the other theoretical results^. 

For low carrier densities, the scattering rate is constant 
and the conductivity acquires a universal constant for low 
temperatures which is twice the quantum of conductance. 
The deviations are quadratic in temperature and carrier 
density. Rather striking is that the hall coefficient is no 



longer scattering rate independent and is proportional to 
the carrier density. 

In the entropy dominated regime (HI) where {|^| << 
fi << T}, the electrical conductivity is linear in tem- 
perature. The scale for linearity is set by the scattering 
rate which in this regime is inversely proportional to the 
chemical potential. In other words, the conductivity has 
a r//i dependence at high T. The hall coefficient again 
is independent of the scattering rate. In this limit of high 
temperatures, for the chemical potential larger than the 
impurity band width, R^ ~ Sn^ec {T/ fi) . 

The hall coefficient is linear with chemical potential at 
low carrier densities crossing over to a ^~^ at large carrier 
densities. The crossover occurs at when the chemical po- 
tential crosses the impurity band width. As such the Rh 
has a peak whose position is sample dependent but very 
weakly dependent on temperature. For ideal graphene, 
where the carrier density would be zero at the node, one 
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FIG. 5: The Hall coefficient is plotted in units of 
a function of chemical potential (/^). The impurity concen- 
tration of unitary scatterers is n^ = 1.7 x 10"*, while the 
temperature of each curve is given in the legend. 



transport but have opposite signs for thermopowcr. In 
particular, thermopower vanishes at the node and falls 
off as [ijT for large temperatures, in contrast to typical 
metals, which it would be a constant equal to the entropy 
per particle of a classical electron gas. 

The calculation of the thermoelectric tensor is tech- 
nically more complicated as the conventional Kubo for- 
mulas need to be generalized to include the effect of 
magnetization^"—. Fortunately the correction due to 
magnetization is cM/T— i^, which in the limit of weak 
magnetic field is proportional to B^ and is neglected to 
leading order in magnetic field. 

The kernels appearing in the calculations for ther- 
moelectric transport are related to those that deter- 
mine electrical conductivity and can be expressed as: 
/3^ = ((/x — e)/e) * o^^. For magnetic field perpendic- 
ular to the graphene sheet the kernels are (see appendix 



would expect the hall coefficient to diverge and change 
from positive to negative as the chemical potential crosses 
zero energy. The fact that in all samples the divergence 
is cut off at some energy scale and the hall coefficient 
is zero at the nodei'^^'^i. For a pure Dirac spectrum, 
the crossover occurs when the chemical potential crosses 
temperature. For unitary scatters, the crossover scale 
is independent of temperature for temperatures smaller 
than the impurity band width. 



VI. THERMOELECTRIC TRANSPORT 
COEFFICIENTS 
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The thermoelectric conductivity tensor i; 



,39. 
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We proceed to analyze the properties in the various 
regimes as before. All expressions are quoted in terms 
/So = ^^£^, wl = 2e\^\vl,/ch and 5*0 = kg/e. 



where the U^f is the correlation function of charge and ^- Thermoelectric transport and Scaling Behavior 



energy current densities. In the presence of an external 
magnetic field, the canonical momentum operators are 
used to define the appropriate current vertices. Since 
two bands touch in graphene, both electrons and holes 
contribute. The two contributions add for off diagonal 



A plot of the numerical evaluation of the thermoelec- 
tric power and Nernst signal are displayed in fig|6]and[7l 
The asymptotic dependences in the regimes identified in 
the previous section are shown in table IIIIl 



J 



r 



At high temperatures, the thermopowcr is related to 
the entropy per unit charge. Since both electron and 
hole states are thermally populated, the net charge is 
an imbalance between the two. At the Dirac point, the 
system is particle hole symmetric and the thermopower 



goes to zeroiSr— . At small carrier densities, the differ- 
ence between positive and negative charge occupations 
is linear in the chemical potential, and the thermopower 
is ex /i/T. This dependence on the chemical potential is 
very similar to the high temperature classical limit. The 
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TABLE III: Thermoelectric coefficients 



Quantity 
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theriTLopower is linear in /i and the coefficient of ^/T is 
a measure of the relaxation time^i^i^. 

Consider the dependence of the thermopower on the 
chemical potential for temperatures smaller than the im- 
purity band width. For small chemical potentials, the 
thermopower grows linearly and reaches a maximum ap- 
proximately at a chemical potential of order of the impu- 
rity band width. For larger values it decreases as T / ^ in 
agreement with the other theoretical rcsults^^i^i^. As 
the temperature is increased and becomes larger than 
the impurity bandwidth, the thermopower qualitatively 
shows the same dependence but the peak now is at a 
chemical potential of order the temperature. In other 
words, as one increases the temperature the position of 
the peak in thermopower as a function of chemical po- 
tential will remain roughly constant until the tempera- 
ture becomes larger than the impurity bandwidth. For 
larger temperatures, the peak will move to larger values 
of chemical potential. 

For chemical potentials larger than the impurity band- 
width, the Nernst signal is proportional to the ther- 
mopower. In this regime, the scattering rate is inversely 
proportional to the chemical potential. Thus for unitary 
scattcrcrs the ratio of the thermopower to the Nernst 
signal is a constant proportional to the applied magnetic 
field. As the chemical potential is lowered and crosses the 
impurity band width the two start to deviate. Within 
this scenario, for a fixed magnetic field, the ratio goes 
to zero as /i/r. The slope of the ratio as a function of 
chemical potential is a direct measure of the scattering 
rate. Within Born approximation, the scattering rate is 
proportional to the density of state which for graphene 
is linear in energy. Thus one would expect a divergent 
Nernst coefficient at the node. This divergence is cutoff 
and the value of Nernst is proportional to t^ . 

In the clean limit, t — > oo, the Nernst diverges as 
the carrier density goes to zero^. The presence of uni- 
tary scatters makes it finite and is proportional to t^ at 
the node. The Nernst signal is still orders of magnitude 
larger than typical metals but the observed values of or- 
der iQ^V/KT implies a scattering rate of order Q.QAeV 
for unitary impurities. The larger the mean free path, 
larger is the Nernst signal and, therefore, one expects a 
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FIG. 6: The thermopower (Sxx) is plotted in units of So as 
a function of chemical potential (n). The impurity concen- 
tration of unitary scatterers is Ui — 1.7 x 10^'*, while the 
temperature for each curve is given in the legend. 
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FIG. 7: The nernst signal (Sxy) is plotted in units of 5o as 
a function of chemical potential (/^). The impurity concen- 
tration of unitary scatterers is n^ = 1.7 x 10~*, while the 
temperature for each curve is given in the legend. 
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TABLE IV: Thermal transport 
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large Nernst signal for large electrical conductivity at the 
node. As we will see later this is a key puzzle in graphene. 



The thermal conductivity can be expressed in terms of 
the thermopower and Nernst as 



VII. 



THERMAL TRANSPORT 



Having discussed the electrical and thermoelectric 
transport in previous sections, we now consider thermal 
transport. With the boundary condition where the elec- 
trical conductivity is zero, the three are related. 
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where K is the energy conductivity tensor which is related 
energy current-current correlation function. 
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where K^^ is the longitudinal thermal Conductivity and 
Kxy is the transverse thermal Conductivity. Given the 
large Nernst signal and anomalous temperature and car- 
rier density dependence of thermopower, thermoelectric 
contribution to thermal conductivity is significant in 
graphene. Analogous to the hall coefficient we define the 
thermal hall coefficient, 



As in the case of thermoelectric conductivities, the ef- 
fect of magnetization is of order B^ and is dropped in 
our analysis. For a magnetic field perpendicular to the 
graphene sheet, the longitudinal and transverse compo- 
nents of thermal conductivity tensor are (see appendix C 
for details)^ 
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The average energy and the specific heat dependence 
on chemical potential, temperature and scattering rate 
are given in section II. We now analyze thermal trans- 
port and all results are quoted in terms of both kq = 



— ^0 



ksTao 7r 



= ^Sobo&ndSo^^. 



Thermal Transport Quantities and Scaling 
Behavior 



The asymptotic behavior of diagonal and off diagonal 
heat transport arc given in table IIVI 



12 






\ 



60 r 



iO--, 



/ 



/' 



.,' '7 






^" 



Temperature 

50 K 

255 K 

400 a: 



Rt 
Rrf 



-0.05 " 




rJ-iWS- 



FIG. 8; The thermal conductivity {K^^) is plotted in units 
of Ko as a function of chemical potential (/i). The impurity 
concentration of unitary scatterers is ra — 1.7 x 10~*, while 
the temperature for each curve is given in the legend. 



Numerical results for thermal transport are displayed 
in figlH] and fig. HI The thermal conductivity depends 
on the correlations of the energy current as well as ther- 
mopower and Nernst. The anomalous behavior of the 
thermoelectric coefficients in graphene has a dramatic 
effect on heat transport. In particular Wiedemann- 
Franz law is not universally obeyed. At low temperature 
thermal conductivity is qualitatively similar to electrical 
conductivity and corrections are of order (Tt)'^. How- 
ever, at the node and for very large carrier densities, 
the Wiedemann- Franz law is obeyed. In the intermedi- 
ate regime, the deviation grows as the temperature ap- 
proaches the impurity bandwidth. In particular a peak 
develops at the node as the temperature is increased and 
becomes of order the impurity bandwidth. 



FIG. 9: The thermal hall coefficient (Rt.h) is plotted in units 
of Ro,T,H ~ I/Bkq as a function of chemical potential (fj,). 
The impurity concentration of unitary scatterers is n^ = 1.7 x 



10" 



while the temperature for each curve is given in the 



legend. In regimes I, II and III it is proportional to the hall 
coefficient with the constant of proportionality being (tq/ko 



'^fc|r/3e^ 



tion. Other relevant parameters arc: (1) band width 
(D = 9.6eV), (2) Lattice constant (a = 1.42^) and fermi 
velocity {vp = 1.0 x 10®m/s). All of the data is taken at 
255K. Given this impurity concentration we can estimate 
the impurity band width to be of order 0.2ey(2000iir). 
Since most measurements are done at low temperatures, 
we are always in a regime where T « \IniT,\. In this 
regime, the scattering rate is very weakly dependent on 
gate voltage up to chemical potentials of order 0.06ey. 



Electrical Conductivity 



VIII. COMPARISON OF EXPERIMENTS DATA 
AND NUMERICAL RESULTS 

A. Unitary Scatterers 

We compare the transport dominated by unitary scat- 
terers with experimental data on graphenei^. Since much 
of the data is obtained as a function of gate voltage and 
temperature, we need to determine the dependence of 
chemical potential on gate voltage. Experimental control 
over carrier concentration is achieved by using a parallel 
plate geometry. For the experimental range of param- 
eters used, we can assume that the capacitance of the 
device is constant which implies that the gate voltage is 
linearly proportional to the number of charge carriers: 

Q = CVg. 

The only fitting parameter is the impurity concentra- 



Thc longitudinal conductivity has a plateau around 
zero gate voltage crossing over to a linear dependence at 
higher gate voltages. This is consistent with a small im- 
purity bandwidth beyond which the scattering rate is in- 
versely proportional to the energy. Since both the slopes 
and the crossover scale is determined by the same pa- 
rameter, the lack agreement is a clear evidence for the 
departure from the unitary scattering dominated scat- 
tering theory. 

The observed asymmetric data is specific to the de- 
vice studied here. The behavior at high carrier den- 
sities is consistent with a small impurity concentration 
of about 1.7 X lO^'*, but the value at the node requires 
an impurity concentration that is an order of magnitude 
smaller. Since the measurements were done using two 
probes, we have subtracted out the contact resistance. 
This is accomplished by realizing that the large gate 
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FIG. 10: Comparison of calculation and experimental data of 
longitudinal conductivity as a function of gate voltage. The 
observed conductivity is much larger than those predicted for 
an impurity concentration of unitary scatterers of 2.2 x 10"* 
and 6.5 x 10"'*. An impurity concentration closer to 1.7 x 10"'* 
is consistent with the data. 
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FIG. 11: The calculated Hall coefficient and experimental 
data is plotted in units of Ro = t§y^ as a function of gate 
voltage. The data is best fit for impurity concentrations of 
1.7 X 10~* for large gate voltages. The overall features are 
sensitive to the impurity concentration of unitary scatterers 
as can be seen by the predicted behavior for n=1.7 x 10~*. 



voltage conductivity is linear which extrapolates to zero 
for graphene. Note that the contact resistance does not 
affect the therniopower and Nernst measurements, but 
does affect the Hall resistance. 



2. Hall resistance 

The fit to the observed hall coefficient is shown in 
figfTTl The hall varies linearly with gate voltage for small 
carrier densities and falls off as l//x^ beyond a scale set by 
the impurity bandwidth. As mentioned perviously, the 
two probe measurement of the Hall data is not reliable, 
and the fact that no agreement with the data is achieved 
beyond qualitative dependence is expected. 



3. Thermopower 
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FIG. 12: The calculated Sxx and experimental data plotted 
in units of So = ^^ as a function of gate voltage. 



The measured thermopower in graphene is a linear 
function of gate voltage for small carrier densities (see 
fig fT2|) . The slope is proportional to r^. By fitting our 
numerical solution to the data, we find that the impu- 
rity concentration of 6.5 x 10~^ (not shown in the figure) 
can account for the slope of the thermopower. We plot 
the dependence of thermopower for an impurity concen- 
tration of 1.7 X 10~^ which agrees with the asymptotic 
dependence at large carrier density. A single parameter 
fit for the entire range is not possible. 



4- Nernst Signal 

The Nernst signal in graphene is shown in fig ll3l It 
is negative at zero gate voltage changing sign for large 
carrier densities. The peak value is large and about 



50fiV/KT. Theoretically it is proportional to t^ and 
is predicted to change sign as a function of chemical po- 
tential and gate voltage. The lack of a single parameter 
fit, as in the case of the thermopower, is obtained for the 
Nernst as well. The value at the node can be fit with a 
smaller value of the impurity concentration, but the cross 
over to the asymptotic behavior is not captured. 

The above analysis suggests that a single impurity con- 
centration of unitary scatterers cannot reproduce all the 
observed data. While no single impurity density fits the 
conductivity, two different values of impurity concentra- 
tion, differing by an order of magnitude, are needed to 
fit the values near the node and the asymptotic behavior 
at large carrier densities respectively. 
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FIG. 13: Experimental data and calculated Sxy plotted in 
units of So = -^- as a function of gate voltage. 



B. Coulomb Scatterers 

It is clear from the analysis in the previous section that 
unitary scatterers fail to accurately reproduce conductiv- 
ity data. Coulomb scatterers have been shown to repro- 
duce conductivity data at large carrier densities. The 
fits to the data with Coulomb scatterers is presented in 
this section. One caveat to note in these fits is that the 
experimental data for Hall and Nernst are outside the 
regime of validity of our theoretical calculations in the 
hydrodynamic limit. In particular the scattering length 
is much longer than the cyclotron frequency. 



1. Electrical Conductivity 

The conductivity data can be reproduced over the en- 
tire range from low to high carrier concentration (see 
fig ll4p . The impurity concentration required is ric = 
1.7 X 10"'^. The finite conductivity at the node is a re- 
sult of the self consistent treatment of the impurity po- 
tential. Even for weak potentials, the induced impurity 
states provide finite conductivity and screening at the 
node. The agreement with data suggests that Coulomb 
and not unitary scatterers are the predominant source of 
scattering in these systems. 



2. Hall Resistance 

The fit to the observed hall coefficient is shown in 
figflSl The hall varies linearly with gate voltage for small 
carrier densities and falls off as l//x^ beyond a scale set by 
the impurity bandwidth. The lack of the agreement is ap- 
parent, but the same drawback of two terminal measure- 
ment precludes any quantitative conclusions from being 
drawn. 



FIG. 14: The conductivity data and calculations for a charge 
impurity concentration of Coulomb scatterers of ric = 1.7 x 
10"'*. Coulomb scatterers provide excellent quantitative and 
qualitative agreement in the entire range of carrier densities 
measured. 
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FIG. 15: Hall resistivity data fit to calculations for a charge 
impurity concentration of tIc = 1.7 x 10^'', and is plotted in 



units of Ra) — 



\B\e' 



as a function of gate voltage 



3. Thermopower 

The overall agreement between Thermopower data 
and the predictions from Coulomb scattering dominated 
transport is significantly better than that for unitary 
scatterers (see figfT6|. The peak's position and slope are 
off by a factor of 0(1) {^ 2). Since thermopower is sen- 
sitive to higher derivatives of the scattering rate with 
respect to energy as compared to electrical conductiv- 
ity, this reflects the difference in the dependence of the 
imaginary part of the self energy of unitary and charge 
scattering potentials (see figlHand figE}. In particular, 
the self energy in the former varies by a factor of '-^ 4 
while the latter changes by a factor of ~ 2 for a change 
in gate voltage from to 0.1 eV. The weaker dependence 
of Coulomb scatterers provides a much better fit to the 
data. 
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FIG. 16: Thermpower data fit to calculations for a charge 
impurity concentration of ric — 1.7 x 10"''. 
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FIG. 17: Nernst signal data fit to calculations for a charge 
impurity concentration of ric = 1.7 x 10~ . 



^. Nernst Signal 

The observed peak in the Nernst signal at the node 
is smaller by an order of magnitude than that expected 
from charge scattcrers (sec fig ITT)) . This result is con- 
sistent with the observation that charge scatterers are 
not sufficient to accurately reproduce thermopower data, 
since both depend on the variation of the scattering rate 
as a function of energy. Furthermore, the crossover to 
the asymptotic behavior is not captured by the Coulomb 
scattering phenomenology. In particular, the positive 
peaks arc not seen in the data. 

A simplifying assumption made in these calculations 
is that the self energy is independent of momentum. In- 
cluding the momentum dependence similar agreement 
with conductivity^! and Hall coefficient^ has been ob- 
tained for finite range scatterers. The inclusion of mo- 
mentum dependence also yields better agreement with 
thermopower data^ while the discrepancy with Nernst 
data cannot be resolved with this generalizatior*^. 



IX. CONCLUSION AND 
ACKNOWLEDGMENTS 

In this comprehensive study of the transport coeffi- 
cients on graphene, we show that the theories assuming 
a single scattering mechanism, either charged Coulomb 
or unitary scattering, are incapable of accounting for all 
the observed data. A number of studies that focus only 
on electrical conductivity emphasize that one or the other 
of these two mechanisms as the source of the minimum 
conductivity and the linear in gate voltage dependence. 
Since the only fitting parameter is the impurity density, 
our calculations show that charged coulomb rather than 
unitary scatterers fit the observed conductivity data. In 
particular if one chooses to fit the minimum conductivity 
using unitary scatterers, the slope at high gate voltages is 
overestimated, while a fit to slope at high gate voltages 
underestimates the minimum conductivity. Given this 
observation one would conclude that unitary scatterers 
are not the dominant source of scattering in graphene. 

To check the validity of the conclusion based on elec- 
trical conductivity, we fit the data for hall resistance, 
thermopower and Nernst all obtained on the same sam- 
ple. This is important because fitting data on differ- 
ent samples require choosing different impurity concen- 
trations which render a quantitative comparison invalid. 
Moreover, the dominant scattering mechanisms in differ- 
ent samples need not be the same. Analytic results in 
asymptotic regimes are provided for all transport coeffi- 
cients, and numerical calculations are used to fit the data 
over a wide range of gate voltage. 

The most striking disagreement is in the Nernst Data 
near the node. Both types of impurities overestimate the 
value by an order of magnitude. The Nernst Data asymp- 
totic values are recovered for the unitary scatterers while 
the Coulomb scatterers over estimate these values by two 
orders of magnitude. In other words in the regime where 
the Coulomb scatterers fit conductivity data, the Nernst 
is significantly overestimated. The inability of a single 
scattering mechanism to get both the behavior near the 
node and at large carrier densities for all observed trans- 
ports coefficients is the principal conclusion of our work. 

Finally, a few comments on the validity of the theoreti- 
cal calculations are merited. We have used the SCBA for- 
malism to compute the transport characteristics. While 
the approximation fails in the vicinity of the node due 
to interference effects, the regime of validity is large as 
seen in the agreement of the density of states determined 
within this approximation and exact numerical methods 
[see [3J and references therein] . A second caveat is near 
the node the effect of the magnetic fields cannot be cap- 
tured within the semiclassical treatment presented. The 
Landau level splitting becomes larger than the Fermi en- 
ergy for a field as low as IT field when doping level is of 
the order 10~^. For higher chemical potentials, the trans- 
port calculations are reliable. However, within a small 
window near the node both the SCBA and the magnetic 
field effects need to accounted for. 
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Appendix A: Electrical conductivity tensor in a 
Homogeneous Magnetic Field 



and derive conductivity tensor in the presence of a fi- 
nite range disorder. The method developed by Kho- 
das and Finkel'stein to calculate the Hall coefficient is 
used^. The conductivity is calculated within linear re- 
sponse. The electrical current density is^ 



In this appendix, we define the electrical current den- 
sity in the presence of a homogeneous magnetic field 



J 









(Al) 



The electrical current-current correlation function is 



K:pir^,rr,r) = < r.( lim ^[(-^^'^, ^ e'l^,) ^ {~^^?^ + e'l^J]^^ in,T)^{rlr')) 






r,~^rf ZTfi f J J J 



> 



(A2) 



In the presence of a magnetic field the Greens function 
is defined as G'(ri,r2,r) = ej:p(^<i>(ri,r2))G'ri-r2(T), 
where G^ -r^ (t) is the gauge invariant Greens function, 
and the exponential factor accounts for the phase ac- 
quired by the particle along the path r ^ — > r f^. 



Choosing the vector potential as A = ^B x r and us- 
ing the form of the phase factor of the exponential of the 
Green's function <I'(ri,r2) = J^^ A{ r )dr, ean lA2l can 
be written as 






naf/ in,rf;T)+ W^'^f (r, ,rf-T)+ W^'^f {r, , r/ ; r) 4- W'^^p^in ,rf]T) 



£j[H^n-^x(^-?|)r(z^., 



+ — X (?t - T]-)fGr^^r,{T)]{Gr,-rA-r)) 



Am 

„2 



Ajrfi 

o2 






e 



2c 

3^ 



2c 



(A3) 



4„,2 ii-^'^r. + ^ X (^ - r})rGr,-rA-rmi^r, + ^ >< (^ " ^f)f {Gr.-r ,{t))] 



(A4) 



Letting R = r i— r f, taking a fourier transform of each we get 
of these correlation functions and labeling K = p + q 
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,lWn 



2 -^ -^ 



r 



To first order in magnetic field (assumed to be in tlie are 
positive z direction) the diagonal and off diagonal terms 



(A6) 



Tree/-?- r, ■r^^ ~e^i\B\ M-^ .. dG^{iWn + in) ~ ~ dG^{iWn)^ 

+ ^^ E (^^-){^ G^i^u.,.)~G^i^u... + .n)^-—} (A8) 



r 



The current-current correlation tensor and impurity self bands is 
averaging Green's function^^ of graphene which has two 

I 



Ul'%{t ^0;^n) ^ ^ Y {(^|,JG^A,^(*^n + *f^)G^^,^(*u;„)} (A9) 

P ,iWri, 

n^f,(^^0;il7) = ^^^J I ^ ^^^fGAA^^{^w^ + ^n)Gc,^{^Wn+^n)GAA,t{^w„) 

P ■,'i'Wn 

+ a E C>J't)^C!AA,ti'^'^-^ + 'i-^)^AA,^{iw„)Gc,jt{iw„) (AlO) 

cpm _f — ' '^ '^ '^ '^ 

G(AA,i3S)(i?,i£) = ^'^^ " ^^'^^^^ ^ (All) 

(ie-E(i£))2-|</,(Z)|2 

Gc(]?,is) = 1'^^^^' ^ (A12) 

(ie-S(i£))2-|</,(:^)|2 
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Performing the sum over frequency and taking the real 
part of the electrical conductivity tensor we get 
I 



^ t 



/ ^ r^'^ l'^' + "^ )lf%JmGAAAe)ImG^AAe + ^) (A13) 



^N..e^0H ^ r ^(,^(,)(s^.) _ 42.a)) „ ^^(^ ^ ^)(5.(M) _ 5.(2..))) (A14) 



p 
Elf) = i?eG^^^^(e + f]){/mG^^^^(e)EeGc,^(e) +i?eG^^.^(£)/™Gc^^(£)} 
J:g''> = i?eG^^^^(e){/mG^^^^(£ + n)ReGc^j^{e + n) + ReG^^^^ie + n)ImGc,^{e + n)} 



(A15) 



r 



Integrating over momentum and taking the limit of an limit is 
ideal Dirac spectrum the diagonal conductivity in the dc 

I 



fT^C' = 



I 



where A and B are defined as e - i?eS(e) and /mE(e). {ReG%B - /mG|g) ^-^^f^ + 2ReGBBlmGBB^^§r^) = 
The hall term can be simplified by letting f2 -> jjjitQj^ 9G^\ 
and using {ReG%g + ImG%g)ImGc = /to(|Gb_bPGc), 



""^v = -^ 2^ j ^^{—Q^i\GAA,'f(^)\ ImGc^-f{e))-nF{e)Im{Gj,j,.^{e) ^| )} (A17) 



r 



Calculating the angular integral first the hall term re- 
duces to 



J 



jjc NsNye^\B\vl f f de dnpis) ~ ^ ~ ~2 dGc .^ ,.^^. 

''''y ^ 2^ / ^/ V^~d (l'^^'4,it(e)l IiTT-Gc-j^ie)) - nF{e)Im{G^^^^{e)-^)} (A18) 

where the integrals of the green's functions are 

~ ~ —IB^ — A^B^+A^ 2AB 

dx{\GAAA£)?ImGcA£)) = ^AB ^ B^ + A^ ' 2AB ''''''^'''' B^ - A^ ^ ^^^^^ 



Integrating by parts, 
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' xy 



-NsN^e^0\vl fdednpie) 1 ^ B^ ^ A^ B^ + A^ 2AB ^ 

'iVTSy 02 , A2 o/iD arctan — ^ —) - 



2CTT 



AB 



TT 9/i ' SAB ^B^ +A^ 2AB 



B^-A^' 3(A2+B2)2 



} (A21) 



Appendix B: Thermoelectric conductivity tensors in field. We follow the same approach as we did for longitu- 
a Homogeneous Magnetic Field 

In this appendix, we define the electrical and energy cur- 
rent density in the presence of a homogeneous magnetic 



dinal conductivity in the previous section. The currents 
are^^ 



J 



,„(r„i) = ]im ^[{-z^^,-e'l^,)^i-z^^^+el'^J]^\n,t)^{r'„t') 



t'~it 



jEAn.t) = lim :^[{-^% - eX")| + {-i'^'^^ + el'^J-^]^\n,t)^ir'^,t') 



r'.^r, 2m 



dt' 



(Bl) 
(B2) 



Letting t -^ it the energy current becomes 



r'-^vi 2m ' ' OT ' * or 



(B3) 



The correlation function that determines the thermoelectric response is 



( hm :^[{i^t - etf, ) - (z^f'^ +eX^)]vI/t(r^,0)vI/(r^,0)) > (B4) 



The greens function G'ri-r2(''") is found using the im- 
purity self averaging technique^^. Since we have terms 
that depend on the derivative with respect to r, we use 
the equations of motion to determine the greens func- 



tions. Using ^G'(ri,r2,±r) = ±S{T)S{ri - r2)+ < 
Tr|r*(ri,T)*t(r2,0) >, eqnlHis 
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nfj(''i,''/;T) = nfj'^(ri,r/;T)+nfj'^(ri,r/;T)+nf'^'^(r„r/;T)+nf'^'''(r,,r/;T) (B5) 

r[(*^n + ^ X (F) - Ft))"G.,.^.,(-r)][(*^., _ 1_ x (F) - Ft))^(<5(r)5(r, - Vf))] 



4to^ 



In ixLomcntum space, with the definition R — p + q , 
we get 






Am?j3 J-^ ^^ 2c '^' ^ 2c 

p - 



nfj'2(^;zl7) = -^ ^ (,^„+,J7)G^(,^„ + ,f7)[(-^ + ^^x^^r(^ + ^^x^y)'^G^(^^„)] 



4m2/3_;^ ^ " ."Kv -" ' '^^ " ' 2c ^^ ^" 2c 

") .iWn 



ei 



Am^(3J-^ Lv ^ . 2c p^ v^ ' 2c 

U^'f{t;^n) = ^ Y. [(-i?-^^x^^rG^(^u;„ + ^f})][(-^-^l-x^^)^H^«„G^(^^«„))] 

P j'ilfn 

- £^ E [(-^-^^x^^rG^(^u;„ + ^f^)][(-l^-^^x^^)^(l)] 

p .zi^Ti 

P ,iuJn 

+ i^ E [(^-*^x^^rG^(^^„)][(]^-^^x^^)/^(l)] 

P ,^''i'n 
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To first order in magnetic field, 



in^ 



nf;,^(V-^0;zf7) = -^ E {((*^n + ^)^| J%(*^n + *")]%(*^"n) + ^| 



G^{iWn + i^) + G-^{iWn) 



} 



nf;,^(^^0;zf}) 



4c/3to 



(iwn + -z-)[v^^^){ — G^[lWn) - G^(iw„ + ^i2) — } 



in. 



+ -Tz^;r E (^"''^ + -T')^~^t,v){ 



dpx 
.dG-^{iWn + in.) 



Acl3m 



P ^iWn 



dpy 



G-^{iWn) ~ G^{iwn + in) 



dpx 

dG^iiWn) 
dpy 



+ 



E {"^^^(*^" + i^)G^{iWn) - G^iiwn + in) + G^iiWn)} 



AcPm J-^ 2 



} 



(B7) 



In tlie thermoelectric correlation function, only the terms conductivity tensor is equal to the electrical conductiv- 

of the form (iwn + ■^) contribute to the real part of ity kernel multiplied by energy divided by the electron 

the thermoelectric conductivity. Performing the sum charge. Thus, the thermoelectric conductivity in the dc 

over frequency and taking the dc limit the thermoelectric limit is 



J 



oDC 



DC 



r'xy 



1 B2_^2 



-NsN^e^\l^\vl fdednp{e) 



2c7r J IT dfj, 

NsNye^ll^lvl f de 



2CTT 



■ / -^nFie){ 



-ie-fi){ 

AB 

3CA2~KB2J2 



SAB^B^ 



} 



B^ + A^ 2AB 

arctan 



AB 



A2 2AB 



B^-A^' 3(^12 + 52)2 



(B8) 



(B9) 



Appendix C: Heat conductivity tensors 

The energy-energy current correlation function is 



K^in^rf;n,Tf) = <T.(lim -^[{-^^^, - et^,)-^ + {-t^^^ + e-l'^J^]^b\r,,n)^irl,Tl)) 



r[^ri 2m 



dn 



dr' 



J J J J 






Using |: G(ri,r2,±T) = ±5{T)5{ri - r2)+ < in r and |^G(ri,r2,T) +S{ri - r2)^S{T) =< 



Tr^*(»'i,T)*U''2,0) > for first order derivatives Tr 



^2 



— Qr'^^ Q^'^^ > for second order derivatives in r. 
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nf;|^(?-'o ?■/; t) = nf ;|^'^(r,, r/; r) + nf '|^'^(r,;, r/; r) + IlfJ'^{n,rf;T) + nf ;|^''^(n, r/; r) (C2) 



[^(-^l^. - ^/)^)[(-*^n - ^ X (^ - ?t))"(*^., + ^ X (^ - ?t))^G.,_..(-r) 



1 

+ i^[(-^^n - ^ X (^ - ^))"(<5(-r)<5(., - rO)][(^^., " ^ >< (^ " ^))^(^%^)] 
+ i^[(-^^n - ^ X (^ - ?t))"(<5(-r)<5(r^. - rO)][(*^., " ^ ^ (^ - ?t))^(5(r)J(r, - r/))] 



In monicntum space, with the definition A = p + q , 
we get 
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(^;»n) = ^ Y. ((*^n + *f^)'G^(^^"n + *f^))[(-l^ + ^x^^r(-^-^x^^)^(G^(*«;„))] 
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To first order in magnetic field 
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By performing the sum over frequency and taking the de divided by the electron charge squared. Energy conduc- 
Hniit of the energy conductivity tensor, the result is just tivity in the dc limit is 
the electrical conductivity kernel multiplied by energy 
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